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We study the binding energies of spin-isospin saturated nuclei with nucleon number 8 < A < 100 
in semiclassical Monte Carlo many-body simulations. The model Hamiltonian consists of, (i) 
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nucleon kinetic energy, (ii) a nucleon-nucleon interaction potential, and (iii) an effective Pauli 
potential which depends on density. The basic ingredients of the nucleon-nucleon potential are, a 
short-range repulsion, and a medium-range attraction. Our results demonstrate that one can always 
expect to obtain the empirical binding energies for a set of nuclei by introducing a proper density 
dependent Pauli potential in terms of a single variable, the nucleon number, A. The present work 

t-h ; 

shows that in the suggested procedure there is a delicate counterbalance of kinetic and potential 
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energetic contributions allowing a good reproduction of the experimental nuclear binding energies. 
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This type of calculations may be of interest in further reproduction of other properties of nuclei 
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such us radii and also exotic nuclei. 
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In our exploratory work we studied the nuclear binding energies for medium mass 
nuclei with nucleon number 8 < A < 44 in semiclassical simulations via Monte Carlo many- 
body techniques. The purpose of that work was to study the role of an effective Pauli 
potential, which is often adopted in semiclassical simulations of many-nucleon systems. It 
was demonstrated that the empirical binding energies for these nuclei can be reproduced 
satisfactorily using the Pauli potential, where the density dependence is parameterized by 
one variable, the Fermi momentum. The agreement with the empirical binding energies 
was excellent in spite of the simplicity of the model. The conclusion of our previous work 
is rather general, it does not depend on the detail of the NN potential provided a short- 
range repulsion and medium-range attraction are included. One can always expect to find 
a proper counter balancing density dependent Pauli potential to reproduce the empirical 
binding energies. Although the model is not based on the fundamental physics of strong 
interaction QCD, it gives a possible guidance for treating complicated many-nucleon systems 
in a simple, practical manner in semiclassical simulations._This may be very helpful to study 

_5[] and neutron halo nuclei 



6]. 



the many-nucleon systems such as the pasta phase J3, [3, 

In this study, we are able to extend our previous work [1] to treat a wider range of nuclei 
with 8 < A < 100 for spin-isospin saturated Z = N (even Z and N) nuclei. We show 
that the density dependence of the Pauli potential can be well parameterized in the whole 
range in terms of one single variable, the nucleon number A. The Fermi momentum is no 
longer a good parameter due to the fact that the average Fermi momenta in nuclei with 
A > 50 saturate to a value ~ 260 MeV/c (e.g., 265 MeV/c for a 208 Pb nucleus) fl. The 
present result generalizes our previous conclusion that the density dependence of the Pauli 
potential is crucial to reproduce the empirical nuclear binding energies. It is important 
to explore heavier systems since stability of this type of procedures must be tested with 
respect to increasing number of nucleons. We have shown in Ref. [l| that in the limit 
of infinite symmetric nuclar matter in these many-body simulations one should obtain the 
corresponding binding energy around —16 MeV. Also in other simulations as in heavy-ion or 
pasta phases one usually has a wide mass distribution of clusters. In addition, we explicitly 
show that it is possible to reproduce the empirical binding energies using different NN 
interaction potentials. Thus, it suggests a simple and pragmatic procedure in modelling a set 
of nuclei calibrated by the empirical binding energies for a given NN interaction potential. 
Then, the evaluation of further experimental observables of nuclei such us radii and also 
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exotic systems such us hypernuclei could be attempted with this type of procedure 

In the present approach, nucleons are treated as classical, structureless particles. The 
model Hamiltonian consists of nucleon kinetic energy, NN (V/vat), Coulomb (Vcw) and 
Pauli {yp au ii) potentials. The Pauli potential simulates nucleon fermionic nature using the 
Gaussian form introduced by Dorso et at [9], but we allow for a density dependence. In this 
study, we use a simplified NN interaction potential keeping only S-wave interactions without 



isospin nor spin dependence [l]. The model Hamiltonian is given by, 

A p2 A 

H = J27r J —+ H l V NN(rij) + V C oui(rij) + V Pau ii(r ij ,p ij )], (1) 

i=l Zm N i=l,j>i 

where pj is the 3-momentum of i-th nucleon and = |r» — r, | {jp^ = |p$ — p 3 -|) the relative 
distance (momentum) of the z-th and j-th nucleons. Explicit expressions for the potentials 
in Eq. ([1]) are as follows. 

• NN interaction potential: 



V NN {r i:j ) = < 



Vcore, for < Tij < a, 

-V , for a < Tij < b, (2) 
0, for a + b < r^. 



The potential consists of a repulsive core of strength Vcore of width a and an attractive well 
of strength Vq and width b. The values used are, Vc or e — 10 MeV, a = 1 fm and 6 = 2 fm. 
For Vo, we use two values, V = 3 MeV and V — 5 MeV. 

• Coulomb potential: 

Vcoulinj) = ^(1/2 +7-0(1/2 + 7}), (3) 

where (tj) is the isospin third-component of i-th (j-th) nucleon (+1/2 for protons, —1/2 
for neutrons), and e the proton electric charge. 

• Pauli potential: 

( r 2 - p 2 - \ 

Vpaunin^pij) = V P exp \-7^2 ~2p\) <W<W;> ( 4 ) 

where 5 TiTj {5 ai(7j ) is the Kronecker's delta for the isospin (spin) third-component. It prevents 
nucleons from occupying the same phase space volume when they have the same quantum 
numbers. (See Ref. |lo| for other approaches.) As demonstrated in Ref. [l], it is crucial to 
allow a density dependence for this Pauli potential if one wants to reproduce the empirical 



binding energies. Thus, we discuss next the density dependence of the Pauli potential before 
presenting results. 

First, we show in Fig. [T] the average Fermi momentum versus the nucleon number A 
obtained by interpolating the values given in Ref. [3|. One can see that the Fermi momenta 
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FIG. 1: Fermi momentum interpolated vs. nucleon number A. 



(Pf) increase as the nucleon number A increases up to A ~ 50. For heavier nuclei with 
A > 50 they saturate to a value of p F ~ 260 MeV/c (e.g., 265 MeV/c for a 208 Pb nucleus) 
Thus, including also heavy nuclei with A > 50, it is more convenient to use the nucleon 
number A to characterize the density dependence of the Pauli potential, although the Fermi 
momentum was used previously [l| for nuclei with 8 < A < 44. In this work, we need 
to extend the parametrization to heavier systems to parameterize the density dependent 
strength V P in the Pauli potential Eq. (j4]). 

For q and p in the Pauli potential, the density dependence can be determined as follows. 
In a nucleus, a typical nucleon sphere radius r may be given by, 



1/3 



(5) 



where p = 2p F /3n 2 is the nucleon density and pp the nucleon Fermi momentum. Then, 
the average inter-nucleon distance 2r may be estimated as (2r/y / 2q , o) — 1 ; where q is "an 



effective range" of the Pauli potential. With the uncertainty principle, go Po — % this leads 
to: 

(9vr) 1/3 ^ 

Po - — = , a s 1/3 Pf- (7) 

For the Fermi momentum pp appearing in Eqs. (jBj) and (j7|), we use the value as shown in 
Fig. [Tjfor the nuclei with 8 < A < 44, while for 48 < A < 100, we use the saturated value, 
p F = 260 MeV/c. 

By performing simulations to reproduce the empirical binding energies for the nuclei 
with 8 < A < 100 for Vo = 3 MeV in Eq. (|2J), we get the A dependence for V P in the 
Pauli potential as shown in Fig. [2] (the blobs). This justifies that the density dependence of 




FIG. 2: Pauli potential strength Vp vs. nucleon number A obtained with V$ = 3 MeV in the NN 
potential Eq. ([2]). Simulations are performed for nuclei with A = 8, 12, 16, 100 (even Z and 
N) spin-isospin saturated nuclei. This applies for all results treated in this study. The solid line 
corresponds to the parameterization given by Eq. (jHJ). 



the Pauli potential may be parameterized well by the nucleon number A. As expected, the 
strength Vp of the Pauli potential increases with A. This behaviour may be analogous to 
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that of the vector potential in Hartree approximation in relativistic mean field models 11]. 
Then, for Vq = 3 MeV, we get the parameterization: 

V P (A) = -25.645 + 2.9596A + 0.0551A 2 (MeV), (8) 

while for Vo = 5 MeV, we get: 

V P (A) = -1088.2 + 140.55A + 0.9809A 2 (MeV). (9) 

Both cases reproduce the empirical binding energies well. Note that the parameterizations 
are approximate and are given as a guidance. Thus, this suggests that, for a given reasonable 
NN interaction potential, we can model a set of nuclei which are calibrated by the empirical 
binding energies. Then, we can use them to study other properties of these nuclei, such 
as charge distribution, proton and neutron density distributions, and proton and neutron 
r.m.s. radii [8]. Although this procedure may be simple, one can regard that all complicated 
many-body effects are condensed into a density dependent effective Pauli potential. The 
present approach is not based on the fundamental theory of strong interaction QCD, but we 
would like to emphasize, a simple, pragmatic aspect for treating complicated, many-nucleon 
systems in semiclassical simulations. 

Now, we are in a position to discuss the results. In the present study, all simulations 
are performed with a fixed temperature T = 1 MeV. In the simulation a nucleus is con- 
structed by initially placing A nucleons uniformly inside a sphere of radius Rq of range 2-3 
fm within a cubic box of volume V = L 3 and impose L >> r^. Then, using the Metropolis 
algorithm [l^ . the ground state configuration is searched by thermal relaxation. The Pauli 
potential should be gradually turned on to avoid instabilities. After this, we sample the 
configurations in order to calculate the statistical average for its binding energy. 

In Fig. [3]we show the binding energy per nucleon (-E/A) as calculated with the model 
with Vq = 3 MeV versus nucleon number A. The central bars show the calculated values of 
nuclei considered in this work and match those in Ref . 13j . The dotted line is a guidance for 



eyes. The statistical uncertainties are shown by error-bars, and they are less than 5 %. The 
empirical values are well reproduced with the density dependent Pauli potential. Thus, for 
this NN interaction potential with Vo = 3 MeV in Eq. (J2J), we have obtained a set of nuclei 
which reproduce the empirical binding energies. Let us note that shell effects are partially 
included through the A dependence in the Pauli potential strength parameterization. Next, 
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FIG. 3: Binding energy per nucleon, —E/A, for the set of nuclei considered in this work. (See also 
caption of Fig. [2J) 

we discuss the dependence on the NN interaction potentials by comparing the results with 
Vo = 3 MeV and Vo = 5 MeV. We show contributions from the kinetic and potential energies 
in Fig. H] for the nuclei with 8 < A < 56. The solid and dashed lines are the results for 
Vo = 3 and Vo = 5 MeV, respectively. For each case, the upper (lower) line corresponds the 
contribution from the kinetic (potential) energy per nucleon. The dotted-line is the sum of 
the two contributions, for both Vo = 3 MeV and V = 5 MeV cases. This suggests that 
the empirical binding energies can be always reproduced by introducing a proper counter 
balancing Pauli potential, once a NN interaction is specified. 

Further, we analyze how the empirical binding energy can be achieved by considering the 
20 Ne nucleus case as an example. In Fig.[5]we show kinetic (K/A), potential (V/A) and total 
(E/A) energies per nucleon for a 20 Ne nucleus versus the Pauli potential strength Vp. The 
empirical value is indicated by the cross. As increasing the Pauli potential strength Vp, the 
kinetic energy contribution increases, while the potential energy contribution stays nearly 
constant. The increase of the kinetic energy contribution originates from the momentum 
dependence in the Pauli potential, due to the modification in the canonical momenta (or 
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FIG. 4: Kinetic and potential energy contributions to the binding energy per nucleon. The solid 
and dashed lines are the results obtained with Vq = 3 and Vo = 5 MeV, respectively. For each 
case, the upper (lower) line corresponds to the contribution from the kinetic (potential) energy per 
nucleon. The dotted line is the sum of the two contributions for both cases, and goes trough the 
empirical values. (See also caption of Fig. [2j) 

the effective masses) of the interacting nucleons via the Pauli potential. Because of this 
positively increasing kinetic energy, the empirical binding energy is finally achieved. This 
example is for a fixed nucleon number nucleus. In order to be able to reproduce the empirical 
binding energies for a set of nuclei, the similar procedure must be repeated for all the nuclei 
in the set. Thus, one can naturally understand why A (or density) dependence in the Pauli 
potential is necessary. 

For the first time we have demostrated that the density dependence of the Pauli poten- 
tial, crucial to reproduce the empirical binding energies in semiclassical simulations, can be 
parametrized in terms of a single variable: the nucleon number A. Such results overcomes 
previous hipothesis describing the density dependence in terms of the Fermi momentum. 
Once the correct density dependence is included through the Pauli potential, the only other 
necessary ingredient is a reasonable nucleon-nucleon interaction. The procedure presented 
is robust, the simulation remaining stable when incresing the number of particles. The 
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FIG. 5: Energy contribution to the 20 Ne total energy. 
The solid (dashed) [dotted] line stands for the kinetic (potential) [total] energy per nucleon, and 
labelled by K/A (V/A) [E/A\. The empirical value is indicated by the cross. 

validity of the parametrization of the Pauli potential in terms of the nucleon number A 
is made manifest by the counterbalance between the growing repulsion coming from the 
kinetic energy and the incresing attraction of the potential energy. This result opens the 
door to study properties of wider sets of nuclei with a correct parametrization of the Pauli 
potential. Although the procedure presented in this work is not based on the fundamental 
theory of strong interaction QCD, we would like to emphasize that, it suggests a simple, 
pragmatic procedure in modelling a set of nuclei calibrated by the empirical binding energies 
for a given NN interaction potential. Then, using this procedure one may think of further 
testing asymmetric systems by studying other properties of nuclei such as radii or other 
exotic nuclei in semiclassical simulations. 
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